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Significant structural, response of a cylindrical 
duct to an internal pure tone sound field resulting from 
an external pure tone source, located approximately on 
the duct centerline, will arise only when coincidence 
occurs between the natural modes of propagation of 
acoustic waves v/ithin the duct and the natural modes of 
structural vibration of the duct itself. 

The coupling mechanism giving rise to such coincidence 
lies within an assumption of small variations of the source 
location from the duct centerline. This result which 
arose from a theoretical analysis based on a solution for 
the velocity potential within a semi-infinite cylinder in 
the presence of a non-axiul incident plane wave and an 
equivalent modal resonator model of the cylinder undergoing 
principally radial vibrations in a lobar axial] v varying 
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pattern was borne out by the experimental results. 
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INTRODUCTION 

The subject of this thesis lies within the realm of 
the interaction of structures and sound. The various 
approachs to sound and structure interaction problems 
in general fall into one of two categories. 

The first which is a relatively new engineering science 
is one based on the theory of randomness in which classical 
solutions are often bypassed when meaningful statistical 
results can be obtained through other approachs. This is 
a particularly useful analytical tool whenever the structur- 
al vibration of interest and/or the sound field of interest 
exhibit random qualties which allow a statistical analysis 
of the problem. 

The second approach and the one within which the subject 
of this thesis rightfully belongs, is a direct analysis 
of a finely structured problem. In the case of this thesis 
the problem can be stated in general terms as follows; 
the response of a cylindrical duct to an internal pure tone 
sound field, where the sound field is the result of an 
incident pure tone plane wave emanating from a source of 
known location. The problem when stated in such a manner 
is deceptively simple sounding, for in actuality the com- 
plete analytical description of the sound field within 
a cylindrical duct undergoing asymmetrical vibrations has 
not been solved in the literature to date. The majority of 
the papers on the subject of sound within a cylindrical 
duct contain the assumj lion that no significant duct 
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vibration exists. In those papers that take duct vibration 
into account the assumption Is invariably made that the 
vibration, is entirely axially symmetric; a condition 
which when examined in light of the light of the forces 
and frequency range involved and other vibration forms 
possible would be very difficult to bring about. Axially 
symmetric vibrations of a cylindrical structure involve 
almost pure stretching energy and are the modes of vibration 
which they are most likely to resist. It must be noted that 
this assumption takes on more validity if the structure 
in question may be termed a membrane ( radius to thickness 
ratios greater than 50). 

A form of vibration, however, for which cylindrical 
structures are uniquely suited is that involving the lobar 
modes. These are the modes of vibration in which the radial 
dispacement varies in the circumferential direction in 
a cos(n©) pattern (n greater than or equal to two). This 
result first given by Rayleigh (8), is obtainable from a 
theory based on inextension of the middle surface Of the 
cylinder wall. The next degree of sophistication to 
Rayleigh’s theory of inextensional vibration is one which 
allows for the presence of axially periodic displacements 
as well as circumferential ones. 

Thus the first purpose of this thesis was to show 
that if significant vibration of a cylindrical duct could 
result from an internal sound field, then that vibration 
must of necessity be an asymm^tirc one involving the nres- 
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ence of lobar inodes and possibly axial inodes as well. The 
second objective of this thesis was to endeavor to describe 
the coupling mechanism between this sound field and the 
cylinder vibration. To accomplish these ends both a 
theoretical analysis and an experimental analysis were 
undertaken. 

The result of these two analyses show that the form 
of vibration resulting from an internal sound field was 
one made up of lobar modes, and that a coupling mechanism 
based on coincidence of the cylinder structural modes of 
vibration and the acoustic modes present in the internal 
sound field does exist and within limits will predict the 
form of vibration if any that will result from a given 



sound field. 
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PROCEDUR E 

The problem that is posed in this thesis may be 
subdivided into a number of distinct components, each 
a dilineated problem within itself and with essence of 
the final results of the thesis lying within the frame- 
work of the interfacial solution matrix of the component 
parts. The component parts are: 

1, A theoretical analysis of sound propagating 
within a semi-infinite rigid duct, where the sound field 
is the result of an external pure tone plane wave source 
offset by a small angle from the duct centerline. This 
is a classical analytical problem in the area of potential 
theory. The desired result of this analysis is a description 
of the sound field within a cylindrical duct at a given 
frequency. The result is desired only in terms of frequency, 
source location, and the relative strengths of the differ- 
ent modes. 



2. A theoretical analysis of the natural modes of 

of vibration of a cylindrical duct of finite 1c: gth. This 
is a topic with a well-established experiment-'] ly v=lid°tcd 
theory. One of these theories with an approximation based 
on different end conditions was used in this nart of 
the analysis. The desired result of this component of 
the problem was the modal form of vibration of a cylinder 
at a given frequency. 

3, An experimental on-fyoi:; of a stool cylinder of 
the approximate dir -n si one ; r - ft. :• n ] c-..o:th, 6 inch - in 



d:i nmeter 
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a pure tone sound field emanating from a loudspeaker locat 
ed a distance of approximately 3 ft. from an open end of 
the cylinder and slightly ( 2° to 3° ) offset from the 
centerline of the duct axial centerline. The desired 
result of this component was a frequency versus signif- 
icant vibration level display for different source 
locations. 

4. Assumptions are the component part which tie 
together the previous three components. The validity 
of the final results lies solely within the validity 
of the assumptions. Basically assumptions played the 
following roles. 

In the' theoretical sound analysis/experimental 
analysis’ interface they took the form of adapting the 



experimental model to fit the requirements of the 
assumptions of the theoretical analysis. 

In the 'theoretical analysis of cylinder vibration/ 
experimented analysis' interface they played the opposite 
role of adapting the theory to fit the experimental model. 

In the interfo.ee botv.cn l.he two portions of the 
theoretical analysis the assumptions tool: the form of 
superposition rather then adaption. Thus the superposition 
of the two separately conceived theories defined a 
coupling mechanism based on coincidence of the acoustic 
and structural modes at a given frequency. 



- 6 - 



is.the interfacing mechanism between the combined theoret- 
ical results and the ensuing experimental results. This 
interface is examined only in the light of the assumed 
variables, 

6. Conclusions and Recommendations make up the 
component part which decides the exact nature of the 
interfacial solutions, their relative strengths and 
weaknesses and finally of deciding what results are 
basic to the problem as a whole , i.e. has the basic 
purpose of the thesis been accomplished. This is done 
in light of both assumptions and the assumed variables. 



THEORETICAL ANALYSIS 



. T heory of Sound Propagation in Cylindrical Ducts 

A considerable amount of literature has been devoted 
to the genera], problem of the propagation of sound in 
ducts. The general form for the solution of the wave 
equation utilizing potential theory as it is presented 
in the following analysis can be found in Morse (l). 

The general solution for the velocity potential 
within a semi-infinite cylindrical duct of zero wall 
thickness, wh en that potential is the result of an 
incident pure tone plane wave emanating from a non-axial 
source, is given by Noble (2). Noble structures the 
problem for solution by a standard Weiner-Hopf technique. 
In the following analysis it will be shown that a more 
direct solution is possible if the external field is 
neglected and the incident wave is assumed to act as 
a massless asymmetric piston at the source end of the 
eylindri c al du c t . 

A. Solution of the Nave Eouation 



Using the 



define a region 



space coordinate system shown in Figure 
within an infinite acoustic medium, 



I 



o < 0 £ a 
o < z < oo. 



Further assume that within that region that the velocity 
potential , 



',0 ( r? y 

1 v N 1 



) 




\ ' 
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Figure I 

Space- Coordinate System 




„q.. 



satisfies the linear wave equation of acoustics, 



V 2 - 



f yjt- 

Co 



T ~ 0 



( 1 ) 



where V" is the Laplacian operator in cylindrical 
coordinates, t is the time coordinate, and C c 
is the velocity of sound in the acoustfc medium* 

Now assume that may he represented hy the 

infinite Fourier transform, 

co 

s \ iUo , ( 2 ) 



-CO 



where W j Cc") is the frequency transform of the 
velocity potential. Similarly let (\L. l LO N ) he 
given hy the inverse transform, 

co 

expfccof] hx . 



-CO 



Then it can he shown that application of do. (2) to 
Eq.(l) will yield the result. 



(3) 



V" 9L * K 2 - 0 



( 4 ) 



where the wavenumber, K , is defined 
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(5) 



It can be shown that the general solution of Eq.(4) 
for waves propagating in the positive direction 

is given by the sum of all modes of oscillation, 

CO CO 






p :=0 



(6) 



* ^A^(o^)cos-p© vBpa(u^^\r\p©|cxp[.t K n iQ 



where 



P‘ C \ 



are positive integers, 



T P ( ^ 



is 



the Bessel function of the first kind, of integer order 



and 



\ 



\ eya 

r i 



is defined b;/ the relation 







( 7 ) 



TTow if th'- surface , 



f 3 <x \ o <• z < co , 

( l ^ 

is assumed to be rigid then it follows that £. 
must satisfy the boundary condition, 




p - a ,o( 



n 

7Z 



JO 



(a) 
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Nov; the application of Eq.(8) to Eq.(6) will yields 



c) QfprjO) 



( 9 ) 



v;here 



*r, 



^pcj 



is the ^ 



TW 



] 



n 



root of Eq.(9). Thus 
p,c| . Nov; if 



, is defined to be, 



has been determined for all 
a cutoff frequency, “) 

r C 0 f & 

• 1 1 

M \ 2 TT 

then Eq.(7) for the mode wavenumber, In pn , may be 
represented by, 



K n ■ K 



1 \ 2 

g) 

\ / 



Vi 



( 10 ) 



( 11 ) 



where, 



i 



CO 

2'ff 



( 12 ) 



It can readily be seen then that for frequencies belo 
the cutoff frequency for a given, mode that the node wave- 
number is imaginary and therefore the mode is highly 
damped and will not propagate. At frequencies greater than 
the cutoff frequency Eq.(ll) yields a real result and the 
mode will oror again, in the positive g- direction with 
the wavenumber specified by Eq.(ll). 
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B. Determination of and . 

Assume that within the region specified by, 

-CD <> -2 r ^ O 

there exists a pure tone plane wave source such that 
the axis of incidence of the plane wave on the the 
point, 

r- o j 2 = 0 

is specified by two quantities B> and (p o ; 

see Figure II. Further let ^ - , represent the frequency 

transform of the unit velocity potential of this incident 
wave in the region, -00$ 2 $ O . Now if it is assumed 
that X-y must vanish in the region, 

a $ r $ oo 

O $ 

then it can be shown that on the surface specified by, 

osru 

2-0 

L must act as an asymmetric massless r is tor 

U 

on the region, 

o « c « o. 

o s 2 « co 

Determination of and may then be 

accomplished in the following manner. If Eq.(6) is multiplied 
successively by Op 0^0 C.OS f © and ijp ( Ypfl ? ) 'S \ A pO 



-1 






Figure II 
Source Position 




and the result integrated over the plane, ~/t ~ ^ , it 

is seen that a non-zero result is obtained only for the 



surface area specified by, 

o s r < o . 

2=0 

Thus the following relation is obtained 3 



„ = (rdr OpC'ipf) • 

s J o v n 






r * 



Ae ^(pe.o.ui) 



eo^ 



pe 



where , 



1 p~C 



(- - 



*s nr \ p > ©■ 

P \zT 



Vz p * ' 



But assuming that the effect of the incident is to 
act as a massless asymmetric piston on the this region 
then it must follow thwl, 

( P, e ,o ,<b) - ^ . (g © ,o ,03) . o< r f a 



and as a result 



A 






<X 



i- (pAr ) p (y P ,r) 



V , 



n 



n 



P 3 'PI 

sir 

A© (g©, 0 ^ 

o 



C6Sp& (lp) 



S\n pa- 
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Now using the notation of Figure II and assuming that 
0-^ = O , with no loss in generality, then the unit incident 
velocity potential is given by, 



— 

+ *>in cb mi© cos 0 



(14) 



• CXp 



-zcot 



and using the well known result. 



co 



-l V C.OS © 



e . 



" 1 m -ern© 

('0 c o 



Tm W 



— ! 
rn=-co 



where vJ yf) is the Bessel function of the first kind 

of integer order m , then Eq,(14) may be rewritten as 






e^K-Z- eo-s 4> 0 - oO-V) 



CO 



J 



m -crnQ 'X / . \ 

(-0 e, 



m*~oo 

Nov taking the Fourier trenui'ovn o: n Eq* ( 15 ) , 

co 



._(*,*>) - \ e °' Vh A 



at 



-CO 

it can be shown that, 

^. ( r - 2 IT 5 (to-co) c 

’ Z CO 



A Z-K- 2 r CO&CD 

s 1 /a* * 



\ ’ / . v m ® “ V / i \ 

) (- z ) C , J m ( Kf - sm ^) 



In r - - f :> 



(15) 



(16) 



v/here 
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^"(ob-co) i s the unit impulse function. Nov: 
assuming that is a small angle such that , 

CCS ~ 1 

s\n <j> b ~ $ 0 

then it can be shovrn that by applying Eq.(l6) to Eq.(13) 
the following results. 



x n 



C 1 
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n 

PI 

2 TC 



y - 2 TT ( C o - IQ 






CO 




irr>G 



COS p© 



le > c-rTe. OmO^kV 



m--co 



s\n p 0 

Examining the integration over \ w ^ first, it 

can be seen that the result is non-zero only for p " “i. 03 



c\© 



and knowing that \T “ ( _ l'j \f then it can he shovrn that. 

*"fy\ ' * Vvi 



A F1 



a 



Z-fT f(on-to') f far ) 0 P (ko^O 






p°l 



O 

ZTT 



. I ( (-0 ? (zf co<, pe 



"COS p 0 



T^i n p 



ho 



Now upon performing the integration X' (*©( ) it must 

follow that (3 pcj is zero for all 



p.<j 



a > and that 



a p , • <£?!?•*) Kftf f T p cr n o 

f P°l o 



CK 



Now define a quantity , such that. 



|3 = k'o. 4> t 



( 21 ) 



Then it can be shown that the remaining integration 

over ^ Af( ") will yield a convergent infinite 

o 0 

series ?.n rowers of p , The result of this 

t 

integration for the four principle modes of interest 
is given below, (the series has been truncated to yield 
accuracy to the fourth significant figure for p < \ ) 
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"'Ct 



A 2l 



» 7T ^ (oh-to’) ( £>) 

r 

- - 5'(co-o.b') jj 



004^ 



(22) 



( 23 ) 



* s , = 2 nr 5*(co-co')[-079i (?4) 
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Now within the region, 

o$? < 0^ 

^ ^ CO 

the sound pressure PC&.t') is related to the velocity 
potential in the following way, 



P(e-dn) - A'-^C’oP) 

O t 



( 26 ) 



and similarly it follows that the frequency transform 
of the pressure within this region is given "by, 



and it therefore follows that, 



Ob 






±. 

2 IT 



ICO 



L 



_ ^ -2 cCvt; 






(27) 



Now upon application of Eq<=(6) to Eq.(27) it rnw.t 
result that, 

i ( Kf'jx-'dv.yu) 



CD 



00 to 



v(f,-tV2 ) , Jp(v l ^y i n C0S P &e 

p=0 Cj = | -Co 



CO O.LO 



( 28 ) 



V/here the of interest are given "by Eq«.’s (22) 

through (25)- 
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A result to "be needed later in tnis analysis is 
evaluated at . This quantity then mus 
the form, 




(, ■- 03 f - • 0044 p^j 

dr, , = CO V ■ & V54 f A ~ ■ <34 4 



"be of 



(29) 

I 

(30) 

(31) 



(32) 
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il. Theory of Vibration of Hollow Cylinders 

A considerable amount of literature exists on the 
theory of vibration of hollow cylinders, particularly 
in the area of calculation techniques for determining the 
natural modes of vibration. Greenspon (3) has compiled 
a synopis of the vast majority of these theories and has 
compared them on the basis of their variance from the 
exact theory and in their ability to predict experimental 
results accurately. It was found that for the duct dimensions 
and range of frequencies considered in the experimental 
portion of this paper that a large number of theories 
gave excellent results. The theory that was chosen was 
that of Arnold and Y/arburton (4)* primarily on the basis 
that it lent itself to a computerized solution. 

A. The Arn old and Y/arburto n Theor y 

The Arnold and Y/arburton theory is confined to those 
modes of cylinder vibration in which bending and stretching 
of the shell predominate. Over the range of frequencies 
and cylinder physical parameters examined in ih i ^ pa ;-ar, 
the bending energy was predominately the oner 'g :invulv<. d 
in the vibrations and the radial direction of vibration 
the predominate direction of displacement. For simplifica- 
tion the rest of this section will deal only with the 
radial modes of vibration duo to their predominance 
in the lower frequencies. Basically in the lower 
range of natural cylinder f: n .el on of rib; : f : Arnold 



nnv ; 



no I vc - 



, . 1 



V'VO'OA f_ 



o:r 



and V/' r bur ton 
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Rayleigh (8), utilizing Timochenko’s (7) strain relations 
and neglecting the trapezoidal shape of the incremental 
surfaces perpendicular to the cylinder axis. Tnis analysis 
allows for the presence of axial nodes as well as circum- 
ferential nodes s whereas Rayleigh's theory does not. It 
is of note however that their theory gives the same 
resonant frequencies as Rayleigh's for the case of no 
axial nodes. 

A difficulty presented itself in that Arnold and 
Y/arburton solved the above problem where the ends of 
the cylinder in question were freely supported. The 
end conditions of the cylinder which was the basis for 
the experimental portion of this thesis were not freely 
supported but rather approximated the condition of a 
free end and an indeterminate end condition. Using the 
coordinate system of Figure II Arnold and Y/arburton assumed 
that the radial displacement was of the following form, 



U(e,f) = L5 0 cos no cos 



r r\n'-g 
” 21 . 



where is the number of quarter wavelengths in the 

axial vibration pattern and *\ is the number of full 
wavelenghths in the circumferential wave pattern. Thus 
it can be seen that this shape function satisfies the 
condition of circularity at the cylinder ends necessary 
for freelv ounrorted end conditions. As stated above, 
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however, the cylinder which was used in the experimental 
portion of this problem was not supported in this manner, 
rather the end nearest the sound source was a free end 
and the other approximated another free end, that end 
being imbedded in sand of a large grain size variety to 
a depth of several inches. An approximation therefore 
had to be made to adapt the Arnold and V/arbur ton frequency 
equation to the cylinder of this problem. The approximation 
was arrived at in the following manner. 

If it is- assumed that the general form of the 
radial displacement in the axial pattern is given by, 

Ufe) C, COSH ) * C 5 (cos V 

1 XX 2U V XX. 

+ C OSH + C , (j m <223 -V S mn ( 34 ) 

■vCuf - SwH 

H ^ 2U iv. / 

then the requirement of of a free end at C> requires, 
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a?- 



o 





(35) 



( 36 ) 



In order for Eq,( 34) to satisfy Eq.(35) and Eq.(36) it 
must follow that (l z and C/, are identically 
zero, or that, 
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deytVrA (■$ S'A •vSn'tH hTil:A 

xi J 11 s xu 
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Nov; a shape function v/hich yields the some resultant 
frequency equation as the Arnold and V/arburton theory 
and which also approximates Eq.(37) for m largo is 
given by s 

Ufe©’) „ Co^rv© - — ^co$rnTr2 - sm (38) 

71 zi~ 



U< 



ST 



It is seen that for large rn that this is closely approx- 
imated by Eqc(37) vrith, 

f? 



C, = " C 



- rnlt£- 



o 



It is easier to visualize this approximation when 
it is realized that for a large number of nodes a 
" free-free " vibration is equivalent to one described 



by Eq 


.(383 


1 with rn an 
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at a given frequency for either set of end conditions, for 

though the effective length has been reduced by a quarter 
of a wave length, 

L ■ L-(l-i') 

so also has the effective number of wave lengths been re- 
duced by the same factor, 

w>' = ai (i- 

A computerized solution of the Arnold and Warburton 
frequency equation is given in Appendix £> , where 

in the output is represented by the letter 

K. 



B. The Equivalent Resonator for a Structural Mode 
Smith and Lyon (5) outline a method hy which a 
structural system vibrating in a single resonant mode may 
be modeled by an equivalent resonator with an equivalent 
mass, damning constant, and spring constant. 

Consider a cylinder of mean radius, (X , length, , 
and density, « Then using the coordinate system 

of Figure II with ©o ~ & , and the mode shape 

defined in Eq.(38) then the instantaneous displacement 
in the radial direction may be given by, 



. \ COS - SNA mV? : 

L 2k 2U J 

Similarily let the instantaneous velocity in the radial 
direction be represented by, 




(40) 




where. 




(42) 





the shape function over the volume of the cylinder be 



equal to the mass of the cylinder, e.g. , 

^ ^ I Lnl = Ho 

r 



( 43) 



then it follows that, 

C Mn 2 \ 2. ‘ 



It can be seen that the kinetic energy, 1 , assoc- 
iated with a given mode of vibration can be represented 
by, 




i- fib v r ,, ( t ) 



(44) 



In the same sense a modal spring constant may be defined 
to be, 



Tf 






1L 

V4i Vl 



1Y\ 



> 6 



(45) 



and a modal resistance, >*»*r 



, such that. 



\ \ ‘ TVv.'l J ) VO T ^ vnr\ ^ 



(46) 



where in a giv^n vibration Tli- is the total power dis- 



sipated internal to the struct 1 .; re,. 
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I II. The Coupling Parameter 

In the preceding section an equivalent modal resona 
tor was defined "by a set of parameters that describe 
the vibration of the structure of interest at any given 
frequency, CO rr , n , 



Next define an equivalent modal force, ^ , 

such that the instantaneous power supplied to the resona- 
tor is identical with the modal component of total power 
delivered to the structure (5). In the present analysis 
where the power delivered to the structure is that due 
to the internal sound field in the cylinder then the 
identity of instantaneous powers is given by, 



Y,„„(f) = (p(ei-) iLf-t) t) As 



(47) 



where the integral implies the the integral over the 
interior surface of the cylindrical duct. Assuming that 
the pressure P (£ ,-fc) can be expressed as the superpo- 
sition of two parts. 



p(e,*A = ♦ FU 0 (.e,*) 



(48) 



then the first, (.£.7*0 f is that pressure which would 

exist in the cylinder in the absence of structural 
motion, i.e., a pressure of the type given by Eq„(29)* 

The second part, , is that due solely to nor -vanish 

ing structural motion. This decoupling of the two parts o 



the sound field is valid if it is assumed that the coupling 
of other non-resonant modes will not occur. This is partic- 
ularly true for the vibration of structures in air, the 
case being examined here, as the density of the structure 
undergoing vibration is usually many times greater than 
the density of air. If it is assumed that Eq.(48) is a 
valid assumption then it follows that the modal force 
may be decoupled in a similar fashion to yield, 
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Now knowing that the sound field of interest is pure tone 
it can be stated that, 



xsl ~ h 



— Z LOfV*V\ t 



w-m, BY, 
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where 



J W* 



is the total luock.d force, 



also follow that. 



- z uWp 



fe> 



- c ^ 



(51) 
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where l 0 is the amplitude of the internal sound field, 
By application of Eq. (5-0 to aq.(tQ) a cenrl i eg ea r-'-otcr 
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X , may be defined to express the ratio between the 

' rnn 

blocked force and the incident sound pressure on the in- 
terior of the cylinder, such that, 



1 
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run 



is given by Eq.(42) and 






v c 



is given 



by the non time dependent part of Eq.(2_9). Thus Ec«,(52) 
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v/h ere for go n e ) \1 i t y ©o 

has beer reinserted. It if. Lb ■ s’: pZ ' to 

show that the coupling parameter, 7 ima * if> non- 
zero only for H f , in which c: e -k; . ( b’j>) 

takes the form, _ ^ 

2 fsuq^ 
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Further it can he ehovm that significant values of 
corresponding to efficient coupling of the sound field 
and the cylinder vibration are achieved only v.hen , 



IV. Theoretical Conclusions 



Eq‘(54) s the coupling parameter, when combined 
with a solution to the natural modes of vibration of 
a cylinder is the basic result of the theoretical portion 
of this thesis. This result is shown graphically in 

A 'll. 

Figure III. The ordinate variable, {\ , is 

a dimensionless frequency parameter defined by, 



A 



2 ~ 
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( f is the frequency in KC) 

where (T is Poisson’s ratio and £ is Young's 
modulus. The abcissa variable is the axial wave length 
factor, A f defined by Eq.(39)« The solid lines 
represent solutions to the Arnold and Warburton frequency 
equation for the cylinder that was the subject of the 
experimental analysis of this thesis in that they 
represent the loci of of modes of vibration having equal 
number of circumferential nodes. The dotted lines 
represent the requirement of Pq«(55) which may be 



combined with Eq.(ll) to yield the relotion, 



K 







(57) 



and are plotted for the ( p , q) modes: 



(4,1) . Thus it must follow 



^ 



(2.1), (3,l),& 

f i cent coupling will 



th-t 



occur only when the (p*l) dotted liner, 
corresponding p— H solid lines* This 
indicated where it occurs on the Figure 



intersect the 
occurence is 
hy a 



o 
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EXP ERIK Eft T.A.I Ai: ALYSIS 



The purpose of the experimental analysis was to 
show experimentally that an internal sound field could 
excite structural resonances in a cylindrical duct and 
that those modes excited would he lohar type modes. 

The form of the experimental model which was 
created to accomplish this purpose depended to a great 
extent on the assumptions made in the theoretical 
analyses. The first section of the experimental analysis 
therefore deals solely with these assumptions and their 
effect upon the experimental model. 

I. Discussion of Theoretical Assumptions 

The following discussion centers only on the 
four major assumptions of the theory which had to 
he met in order for any valid comparison to he made 
between the theoretical and experimental results. 

A. "Sound field restricted to interior of the 
cylindrical duct. . , " 

Figure IV shows a scale plan of the experimental 
model. The sneaker was surrounded completely by two 
layers of three inch fiberglass blanketing on the verti- 
cal sides and three layers on the top and bottom. The 
whole was hold together by thin wooden studding along the 



outside edges to 
aries were prese 
of the interior 
then freely 



insure that no significant solid bound- 
nted to the sound field with the exception 
walls of the duct. This cage-like structure 
susv -nded from a steel ring a) choree 1 in 



was 



sue. 






Figure T V 
Experiment 0 .! Model 




V C 



the ceiling. The cylinder mouth was exposed to the sound 
field by cutting a hole of the exact size of the cylinder 
outside diameter in the bottom of the speaker enclosure 
and then arranging the two so that the end of the cylin- 
der was flush with the inner layer of insulation and 
snugly fitted by it as shown in Figure IV. It was found 
however that the sound dampening characteristics of the 
insulation were not of an efficient nature for fie_quencies 
below 1900 Hz. Below this frequency there existed a sig- 
nificant external sound field. A cylinder vibration mode 
corresponding to the lowest natural lobar mode of the 
cylinder was excited at a frequency of 720 Hz; however 
it was discovered that significant changes in the accel- 
eration levels being recorded could be achieved simply by 
the movement of persons or objects in close proximity to 
the cylinder. This same phenomena was not present at 
frequencies above 1900 Hz. 

B, " Duct acoustically semi-infinite in length... :l 

This assumption led to the placement of tufts of 



insulation matc-'ial to a depth of six inches in the 
bottom of the cylinder. Based on the experience wi ti- 
the externa] sound field described above it was then as- 
sumed that at frequencies obove 1900 Hz that no signifi- 
cant reflection of incident waves occurred at the bottom 
termination of the cylinder, thus allowing only waves 
propagating in the positive Z direction. 
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C. "Pure tone plane wave source..." 

The input to the driver element of the speaker was 
an amplification of the output of a pure tone oscillator. 
The driver element to duct rnouth distance used in the 
experimental data to be presented was approximately 
three feet corresponding to a value of, (see Figure V) 

~ 2 c ( 

thus satisfying the requirement of K ^ l necessary 
for the assumption of a plane wave source. 

D. "Source offset from centerline of duct by a 

small angle. . . , and ^ 1 . . . . " 

Due to the physical, characteristics of the model 
the offsetting of the source location from the duct 
axial centerline was accomplished in a reverse manner 
by inclining the cylinder to a given angle and then 
readjusting the whole such that the center of the top 
of the duct vras directly beneath the the source. Teasure- 
rnents of the inclined angle were accomplished with a 
plumb bob device. The angular offsets used, wore 2° and 
3° thus resulting in a value for ^ of approxi- 
mately 0.3 i (See Figure VI), thus satisfying the 
requirement of ^ < I 

E. "Cylinder end conditions of 'Free ~ Indeterminate' 

in experimental model..." 

The upper end of the duct was truly a free end in 
that it was in no way sunnorted. by the insulating enclos- 
ure. The bottom end was imbedded in sand to a depth of 



Figure V 

Variation of 
with Speaker Location 
and Frequency 
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Figure VI 

Variation of ft with Speaker Location and Frequency 




FREQUENCY ( KC) 
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several inches. The sand could probably be better termed 
small grained gravel in that it supplied sufficient 
holding force to support small, offset angles when the 
cylinder was inclined but could not have been of much 
support when trying to oppose high frequency oscillations 
on the order of microinches. Thus the end condition for 
this end must be termed indeterminate. 

A table of the properties of the cylinder used in 
tne experiment is given in Table I. 
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Property 
Length 
Mean Radius 
V/all Thickness 
Density 

Young’s Modulus 
Poisson’s Ratio 
Inside Radius 

Commercial Type: 



Table I 

Cylinder Properties 

S ymbo l 

L 

a 

t 

E 

a 



Quantit y 

60 in, 

3.1725 in, 

0.280 in. 

0,283 lbs/in 5 
29. 6*10 6 lbs/in 2 
0.290 
3.0325 in. 



6” ASA SCHEDULE 40/ STD BLACK /A STM 53 



II. Experimental Procedure 



Figure VII shows a schematic representation of the 
instrumentation of the experimental model. For the 
data presented in this pacer the following instrument 
checks were first accomplished with a vacuum tube 
voltmeter. 

1. The output of the oscillator was adjusted to 
a set level, measurement taken at point (2)„ 

2. The combined output of the oscillator and 
amplifier was then adjusted to a set level, measurement 
at point ( 3) . 

3. The voltmeter was then attached so as to 
measure the output of the accelerometer preamplifier 
(gain of unity) which was also an input to the oscillo- 
scope, measurement at point (4). 

The actual measurements of acceleration levels 
was then accomplished in the following manner. The 
accelerometer which used ° magnet mount for attachment 
to the cylinder wall was pi need on the surface of the 
cylinder. The frequency was then v ried w. : th the o~c : Hat; r 
until a significant level of acceleration v/>s obm.ei vcd. 

At this point the output of the nmnlifj or was r echo eke 3 
to insure that it had roained the seme. 

Acceleration readings were then Lap on in millivolts 
(rms) at one inch intervals in the axial direction and 



at intervals o 






i n 



I i i y* 
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Figure VII 

Schematic of Exoeriment Instrumentation 
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corresponding to sixteen equally spaced circumferential 
stations, see Figure VIII. This was accomplished by a 
grid network of tape laid out on the surface of the 
cylinder and marked at the proper intervals. The 
sensitivity of the accelerometer used was a constant 
51 mv/g up to 4000 Hz. Above 4000 Hz. the magnet rnountin 
device caused the sensitivity to become inaccurate and 
as a result the experiment was limited to a high freq- 
uency limit of 4000 Hz. Frequencies v/ere determined by 
the oscilloscope display. 
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Figure VIII 

Identification of 
Circumferential Stations 
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III. Analysis of Data 

As previously stated rras acceleration levels were 
measured with an accelerometer along a tape grid which 
remained fixed throughout the series of experiments 
which resulted in the following data. Six different modes 
of significant acceleration levels were observed. These 
six modes are given -in Table II in terms of the assumed 
experimental variables. 

Tables III through VIII of Appendix A are the tabu- 
lation of the original rrns acceleration levels recorded 
in millivolts (mv). The information presented at the top 
of each of those tables was arrived at after two steps 
of data reduction. These two steps are explained in 
detail in the following section. 

A. Data Re duction Steps (l) a nd ( 2) 

Step (l) Det erm in at i on of Basi c Vibration pat tern 

It was assumed that the circumferential vibration 
pattern was of the form. 

co~> n (© - ©o') 

and sirailarily that the. axial vibv:t? or. v-,ttern v >. r of 
the form,. 

CD'S *35* 

XU 

or equivalently, 

Cos ^ 



n ^ 

/*• l-c 






IL 

4 



dn, 



where /\ 



is 



* 



). 
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Table II 
Observed Modes 



Frequency 

f 



Hz 

1950 

2500 

2300 

5070 

5100 

5850 



Source Location 



0 © 

_ o o 

( )° Sta.# 

5 ° 8.5 

2 ° 5.5 

5 ° 8.5 

2 ° 5.5 

5 ° 8*5 

2 ° 5.5 



HC- 



At first glance it is obvious that there is e per- 
iodic pattern to the data in Tables III - VIII* It is 
an easy though tedious exercise to arrange this data 
in alternating positive and negative value groups. If 
this data so arranged is then transferred to a plot 
of amplitude versus position the basic cosine pattern 
is readily apparent. Define the following quantities, 

A nn : Maximum acceleration (rms) in the total vibra- 

la 3.X 

- tion pattern. 

A : Maximum acceleration (rms) recorded at an axial 

a, max v ' 

station. 

A : Acceleration (rms) recorded at an axial station, 

a v ' 



A 



Cj.max 



A 



c 



Maximum acceleration (rms) recorded at a circum- 
ferential station. 

Acceleration (rms) recorded at a, circumferential 
station. 



Figures IX through XV then represent plots 
acceleration level versus position of mens 
plot represents the data for one of the si,. 



of normalized 
uremen*! . Each 
("odcr given 



in Table II. The position of nodes and ant vnodes can then 
be determined by simple examination of each individual 
plot. The actual determination of the nodes was performed 
on plots of too large a scale to be presented here. 

Table IX of Appendix A is a tabulation of the result of 
this determination. 
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Due to physical obstructions it was not possible 
to take axial measurements over the entire length of the 
cylinder. Thus TYl and as a result 7v > we re assumed 
to be given by the following relation, 

rn = Ilk Lin) 

average node separation (in) 

The value of A was easily determined for each mode as 
it must of necessity be the number of nodes in the circuifl' 
ferential pattern. 



/ 
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Figure IX 

Legend for Figures X through XV * 
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Stop (2) Determination of A fflax , \ my , U max 



The first calculation necessary for this step in 
the data reduction vras to determine the maximum value 
of rms acceleration (mv) for each mode. If it is 
assumed that the vibration pattern is given by, 

A (*,e) = *?) 



then the axial measurements taken at ( © G ) must 

v • o'o 

be of the form, 

A A (2-,(e-ea) - *?) 

and similarity the circumferential measurements taken 

at Z must be of the form, 
o 

A c (Zo 1 ©) = Ac.wio eot, n (e-Oo") cos 4 4 ) 



Now if Z is an antinode then v:c must have that 

o 



A 



max. 



A 



c s na 



Similarity if (G - O ) i r. 
' ' 00 



15 node t 1 ' 



A ~ A 
max a s . max 

Nov; if neither of these conditions has occurred then the simplest 

fanner in which to proceed is as follows. Make art enlarged 

plot of the recorded date nearest the intersection of the axial 

and circumferential mea/'w-'i ; g lines, i . a . ; a t the point d: fined 

by (© ~ © ) . Z r 

^ ' o ' o • o 



Define two angles* 



-pa- 



°^C * \ 

O ^ < £ \ 



Me' 

where is the angular displaccrcnt in degrees of the 

iferential 
axial 



4. ( circumferential \ 
neareS0 \ axial ) 



node from the intersection point 
of the two measuring loci.. In a similar fashion define two 



( <A C ) 0 ) 

amplitudes y > 

K ) 



to be the acceleration level of the 



Circumferential position.) . , , . , .. . j. m . 

Axial Position \ at the ”>tersec cion Point. Then 



max can be closely approximated, by, 

fu 



k 



■v 






mfW 



l *Sm <> tn 



For three of the modes observed it was possible to 
take A directly from the data as an ant inode occurred 
at the intersection point. The other three mode s , though , 
required the use of the above technique to dete-'-mine- A . 

i.i cLj\. 

The following three figures show IM s deter 'in a ''ion 
graph i cal ly. 
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Figure XVI 

Legend for Figures XVI Ithroi gh XIX., 
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ANALYSIS of results 

The first point upon which a;n analysis should be 
made is to . what degree did the experimental analysis 
results prove or disprove the two fundamental purposes 
of the thesis, in light of the assumed variables frequency 
and source location. The two basic purposes of the thesis 
were to show that if significant vibration levels could 
be excited in a cylindrical duct by the presence of an 
internal sound field that these vibrations must be of 
the lobar type, and that if these vibration modes exist 
then there exists a coupling mechanism which will predict 
their occurrence. 

I. Existence of Lobar Modes 

The experimental analysis showed that for the frequen- 
cy range examined that the only modes present with signif- 
icant levels of vibration were in fact lobar modes with 
axially varying vibration patterns. 

II. Validity of t h e Co u pling Mechanism 

The coupling mechanism as described in the theoretical 
conclusions may be examined o.i a number of points. The 
first and obvious cotup°ricon :1s acconpl 1 ishod bp plotting 
the proper information from the reduced d at' 1 portion of 
Tables III through VIII on Figure III of the theoretical 
conclusions. This has been done and the points arc 
indicated by a ^ . In each case the points are 

removed from the uredicted value only by the variance 
in A i except for the mode observed at 1950 Hz which 
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deserves special attention and v;ill be discussed fully 
in the Conclusions. 

The second comparison is the location of the nodes 
and antinodes of each mode, for the theory predicts 
their location also as a function of the source location. 
The theory predicts that the circumferential vibration 
pattern must be of the form, 

cos n (© - Go') 



and thus nodes must exist at, 



0 - ©C " V 



■Sir 

zn 



S = 0 3 j ' ) 



and antinodes at, 



©= e 



o 



-V 

zn 



S'- 0,2, "-.W 



Similarly the theory predicts that the axial vibration 
pattern must be of the form, 






and thus nodes must exist at, 



© ( »•*) s — (zs -t) 

m s y 












and antinodes 



2<>1 - C 2 ^ 1 ) s= 



4 • 



m 

Thus it is seen that a shift in 0 O must result in a 
shift in the circumferential nodal pattern. Figure 20 
is a striking exarnnl e of this oh or on on -a and is a riot 



of the normalized amplitude versus circumferential station- 
spacing for the two modes observed at 2300 Hz. The points 
are data points and the curves are not fairings through 
these points but the actual cosine variation predicted 
by the theory. They are faired in each case only to their 
amplitude. The amplitudes of both patterns have been 
normalized by (tf Q = 3°K 

The other modes have been compared in a like manner 
and a result of this comparison is presented in Tables X 
through XV of Appendix A. 




o 
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The second comparison that can he made in terms of 
the validity of the coupling parameter is wh at effect 
the variation of source location had on the amplitude 
of vibration. This can he attacked in the following 
way. In order to make any type of comparison utilizing 
the magnitude of the coupling parameter some assumptions 
must first he made concerning the modal impedance of 
the cylinder undergoing vibration. The subject of 
the variation in modal impedance with mode form, frequency 
and form of excitation is a lengthy subject and cannot be 
be gone into at great detail in this thesis if for no 
other reason than it was an item over which no control 
was maintained and no means were available for determining 
its effect upon the levels of vibration observed. This 
is primarily due to the real rrart of the radiation- 
resistance a factor strongly couoled with the P , term 
of Eq.(48) which was completely neglected in the theory 



portion of the thesis. However it does not seem too 

presumptuous to assume that for the two p airs of moles 

that were observed at the s frequencies, that the 

impedance remained cor si; mt, i « e. , it in co a L-a. t 

with respect to amplitude of vibration. Thus a comparison 

can be made in the following manner. Assume that Die 

modal impedance Z is proportional to the ratio 

inn 

of the absolute value of the coupling parameter to the 
modal rms imxirrun velocity, (6)? or. 



1 ^ \ 



I X rn \ / r-v GOVOi m-TrOv) t 

V 

Vvav\ 



Now if it assumed that for all tne modes except the 
one observed at 1950 Hz, that the requirement, 

= ™ir 

2 L 

has been met then Eq.(54) becomes, 

qfT ^ C ^ <0 Onn C >< 

where the values of are given by Eq.'s(^O) 

through (32). Thus remembering that, 

f> = 

that, 

\* mx \ - C,(Ka) a 





p 

'-l 









C, (Kaf ' (.«3 






1VC. 

'ia 1 



The value of these relations along with if. other 
parameters needed to make the proper comparisons 
are given in Table XVI. The value for |"X mv ,|for xhe 
mode observed at 1950 Hz was calculated in the follow:- ng 
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Table XVI 



* 


n 


y ft 

% 




(,KcO 


lin 

_c, 


Vrrw\ 


Hz 




radians 
*10" 4 




20.6 


•no" 4 

'ret'zxr- v wn-<«>r»ia 

0.86 


in/sec 
*10" 3 


1950 


2 


27.5 


.486 


1.286 


2™ 


2 


12.1 


. 486 


33.9 


19.00 


3.860 


2300 


2 


27.5 


.486 


33.9 


43.00 


6.610 


3070 


Is 


42.2 

*10~ 2 


.434 


346.0 


5. Cl 


0.628 


3100 


3 


144 .0 
-2 

*10 ^ 


.434 


361.0 


17.85 


1.325 


3850 


4 


147.0 

*10'* 


.399 




4590.0 


17.00 


3.130 
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manner, Eq.(ll) was utilized to determine the mode 

wavenumber for this mode. It is then seen that the 

result of Eq.(ll) is imaginary and therefore it is 

a damped mode. It can then he shown that the following 

- ( K I ) 

approximation is accurate for, e ' v pq J '-~'0, 



L 




As stated before the most valid comparison 
based on these parameters will be the comparisons 
of modes at the same frequency and mode shape. 
Thus examining the two modes at 2300 Hz. first it 



can be 


seen 


that two 


comp 


arisons might be made. 
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- ( 
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2) 


Ex 
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. 650 
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45.00*10~ 4 /6. 610^3 0” 5 ) 


0 = 2° 
y 0 
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* ( 


19.00 • 10“ 4 /3 . 86* lO" 3 ) 
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This same comparison may be also made for the two 
virtually ' identical modes observed at 3070 and 3100 Hz. 

1) Solving for the velocity ratio; 

Predicted : 3.36 (17.85/ 5.01) 

Actual : 2.11 (1.325/ 0.628) 

2) Checking constancy 

$ 0 s. 3° : 1.545 = ( 17 « 85*10”^ /l, 325*10"^) 

0 = 2° : .798 = ( 5.01*10“ 4 /0.628*10“ 5 ) 



Finally the other two modes not 3 ^et examined may b 

\7r,J 



checked also for constancy of 



C, 



]/v 



1950 H z 

3 C 



0 = 



o 



>067 = ( 0.86*10"Vl.286*10“ 5 ) 



3850 Hz 



to' ~ 2 



543 - (17.00-10" 4 /3«130-10" 5 ) 



CONCLUSIONS 



The analysis of results has shown that the two basic 
purposes of this thesis were accomplished. Lobar type 
structural vibration modes were excited in the cylindrical 
duct by a pure tone sound field and the theoretically 
derived coupling mechanism within limits predicted the 
frequency at which these modes would appear, the form 
they would take, the location of nodes a.nd antinodes, 
and to a degree even predicted their relative amplitudes 
of vibration. 

The prediction capability of the coupling mechanism 
however was weak in the area of predicting the axial 
vibration pattern that would be present at a resonant 
response frequency. The two obvious answers to this 
problem are at once evident. 

(i) The experimentally observed structural modes 
were in actuality non-resonant modes.. 

(ii) The predicted theoretical modes were incorrect. 
For statement (i) to be true theft it would s°om that the 
resonant modes of ihe s- ... 3 oh or type should also have 



been observed, experimentally. This did not occur. There- 
fore state;. Slit (ii) seems to be the answer to the problem. 
This is even easier to believe when it is realized that 



the assumptions made in adapting Arnold and Tarburton ' s 
frequency equation to the experimental model were undoubt- 
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of the observed vibration patterns fit the theoretically 
assumed shape very nicely, but this does not imply that 
such an assumed shape should also give the same frequency 
result. A solution to the natural modes of vibration 
of a "Free-Free" cylinder with all end conditions 
properly satisfied would undoubtedly have a better fit 
to the data of this thesis than the solution used. 

As was stated in the analysis of results little can 
be said about the capability of the coupling mechanism 
to predict relative strengths of vibration due to the 
unknown properties of the modal impedance. However 
the comparisons that were made in that section of the 
analysis were far better than would be expected consid- 
ering the margin of slack that must exist to some degree 
in both the mechanism and the observed data and the 
powers that many of these quantities were raised to 
in order to make comparisons. Conclusions, other than 
vibration level increases with increase in source offset 



angle should not be based on the 'oresont data. However 
it must be noted that this was not an initial desired 
result of the thesis. 

Perhaps the most interesting result of the experiment 
analysis was the single non-predictod node observed at 
1950 Hz. This mode would seem to be a resonant cylinder 
mode if it is assumed that the modes at 2^00 Hz were 
resonant modes. Also it is not nee ess rily correct to 
term it a non-pred icted mode for the theory docs predict 
a finite degree of cou cling for nodes not too far b-'low 
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the cutoff frequency of interest. Another possibility 
for explanation presents itself if the mode observed at 
3070 Hz is closely analyzed. For the purpose of analysis 
this mode was assumed to be identical with the 3100 Hz 
mode but in actuality its axial vibration pattern varied 
slightly from this sister mode and this is why in the 
data tables it is attributed with a value of m- 16. 5» 
a not too plausible value but the one that best fits 
it nevertheless o Thus it is possible that this is in 
actuality the same type of mode as that observed at 
1950 Hz. It must be noted that the 1950 Hz mode 
best fit the m -11 vibration pattern which is the 
next level below the m - 1 3 pattern of the 2300 Hz modes. 
Thus the interesting phenomena that would seem to present 
itself, and does not seem too unreasonable, is that the 
cutoff frequency is not necessarily a constant value 
but rather might be of a fluctuating or finite bandwidth 
character, allowing close frequencies of vibration on 
either side of it to sway it from its path so to speak. 
This phenomena is definitely an item that deserves fur i. her 
study, both theoretically and experimentally. 



REG OH I -1 EI'-I D A T J 0 E S 

The logical first recommendation would be to cont- 
inue the experimental analysis over a wider range of 
frequencies and with cylinders of varying dimensions 
and mechanical properties. This should be done in 
conjunction with a better solution to the natural 
frequencies of cylinder vibration, the best approach 
probably being to adapt the experimental model to 
a specific theory, 

A second recommendation which should offer interesting 
results would be to excite the cylinder with noise rather 
than pure tone sound. 

A third, recommendation already stated in the con- 
clusions would be to investigate the cause end occurrence 
of modes of the type observed at 1950 Hz, 
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Aouendix A 



Table III through Table V III 

Tabulation of original data in millivolts (mv) and sumrr.ar3 r 
of mode charaeteritics. 

Table I X 

Tabulation of observed nodes and antinodes. 

Table X through Table XV 

Comparison in tabular form of location of observed 
and predicted nodes and antinodes * 
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Table IV 
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Table VI 
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Table VII 
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15.6 


19.0 


3.9 


22.6 


2.5 


26.0 


6. 5 


29.7 


5.0 


33.0 


9.2 


36.5 


7 o 


40,0 


11.8 


43.2 


47.0 


14.5 


50.1 


13.4 






Table X 



f = 1950 Hz 


n = 11 


0^ at 3ta.# 0.5 


n = 2 


Nodes 


Anti . ides 


Circum, 


Axial 


C i r cum . 


Axial 


Sta„# 


Z-inches 


Sta. fr 


Z-inches 


Actual 


Theory 


Actual 


Theory 


Actual 


Theory 


Actual 


Theory 


2.7 


2.5 


— 


2.73 


0 . 5 




— 


8.17~ 


7.4 


6 . 5 


14.2 


13.61 


4.5 


4,5 


19.0 


19.10 


11.0 


10.5 


25.0 


24.50 


9.0 


0 

00 


30.0 


30.00 


15.1 


14.5 


35.2 


37.50 


13.0 


12.5 


41 . 0 


40.90 






45.0 


46.40 






50.0 


51 . So " 






— — 


57.30 











































































































Table XI 



f = 2300 IIz 


m = 13 


at Sta .# 3.5 


n = 2 


Holes 


Antinodes 


Circum. 


Axial 


Circum. 


Axial 


Sta. jf 


Z-inches 


Sta. 7 j 


Z-inches 


Actual 


Theory 


Actual 


Theory 


Actual 


Theory 


Actual 


Theory 


1.6 


1.5 


— — - 


2.31 


3.6 


3* 5 


6.92 


5.6 


5.5 


— -™ 


11.54 


7.6 


7.5 


16*0 


16.15 


9.6 


9.5 


20.6 


20.75 


— 


11.5 


25.5 


25.38 


13.6 


13.5 


29.8 


30.00 


16.0 


15.5 


34.0 


34.62 






33.9 


39.20 








43.84 






— 


48.50 






— 


53.07 








57.70 
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Table XII 



f = 2300 IIz 


m - 13 ( 


at Sta. // 8.5 


n = 2 


LTOxIea 


Antinodes 


C i r cun . 


Axial 


Circurn. 


Axial 


Sta, // 


Z-inches 


Sta. ff 


Z-inches 


Actual 


Theory 


Actual 


Theory 


Actual 


Theory 


Actual 


Theory 


2.3 


2.5 


— 


2.31 


C. 3 


0.5 





”6792 


6.3 


6*5 


11.5 


11.54 


4.3 


4.5 


17.0 


16.15 


10.3 


10.5 


20.8 


20.75 


8.3 


8.5 


26.2 


25.38 


• 

1 — 1 


14.5 


29.8 


30.00 


12.3 


12.5 


34.5 


34.62 






39.6 


39.20 






44.0 


43.84 






48.4 


48.50 









53.07 






— 


57.70 



























































































Table XIII 



£ - 3070 


m =16.5 


0^ at Sta«# 3.5 


n ~ 3 


„i OuOS 


Antinodeo 


Circum , 


Axial 


Circum. 


Axial 


Sta*// 


Z« inches 


Staef' 


Z-inches 


Actual 


Theory 


Actual 


Theory 


Actual 


Theory 


Actual 


Theory 


”2.1 


“2716 


— 


1,82 


6,0 


0.83 




5. 46 


4.8 


4.83 





9.10 


3.4 


3.50 


— 


12.70 


7.3 


7,50 


15.8 


16.25 


5.9 


6.17 


20.5 


20.00 


10.3 


10.17 


23.6 


23.60 


8.9 


. 8.84 




27.3 


12.8 


12.84 


51.5 


30.90 


11.5 


11,51 


34.1 


34 . 5 


15.1 


15.51 


38.6 


38.30 


14.1 


14.18 


43.0 


41.8 






46.4 


45.4 






— 


49.00 








52.7 









56.4 








60.00 



























































Table XI7 



-f ^ 3100 Hz 


m = 17 


e o at Sta.# 8.5 


n - 3 


;;o:or: 


Anti nodes 


Circum. 


Axial 


Circum. 


Axial 


Sta. # 


Z-inches 


Sta. // 


Z-inches 


Actual 


Theory 


Actual 


Theory 


Actual. 


Theory 


Actual 


Theory 


1.25 


1.84 


— 


1.76 


15 . 6 


0.51 


— 


5.30 


3.90 


4.51 


— 


8.80 


2.5 


3.18 


— 


12.35 


6.50 


7.18 


15.8 


15.85 


5.0 


5.85 


19.0 


19.40 


9.20 


9.83 


22.6 


22.90 


7.9 


8.50 


26.0 


26.50 


11.80 


12.50 


29.7 


30.00 




11.19 


33.0 


33.50 


14.50 


15.17 


36.5 


37.10 


13.4 


13.86 


40.0 


40.60 






43.2 


44.3 0 






47.0 


47 . 60 






50.1 


51.00 






— 


54.70 
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Tahle XV 



£ - 3850 Hz 


m - 4 


at Sta*# 3c- 5 


n - 4 


7o7es 


Antinodes 


Circum. 


Axial 


Circum. 


Axial 


Sta *# 


Z-inches 


Sta, # 


Z-inches 


Actual 


Theory 


Actual 


Theory 


Actual 


Theory 


Actual 


Theory 


L 

Ij * 

O 

i 


0.5 


— 


7.50 


TIT 


~r/5 


21.0 


”*22750 


2,2 


2.5 


37.2 


37.50 


3.2 


3.5 





52.50 


4.2 


4.5 






5.2 


5.5 






6 * 2 


6. 5 






7.5 


7.5 






8.2 


8.5 






9.3 


9.5 






10.2 


10.5 






11.5 


lf\ 

( — 1 
i — l 






12, 2 


LT\ 

cC 

i — 1 






13.3 


13.5 






14.2 


14.5 






15.5 


15.5 
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Apoendi x B 

Computerized Solution of Arnold and V/arburton( 4) Frequency 
Equation for a Thin-balled Cylinder Freely Supported at 



the Ends. 



K' A | > ' 



n A- ] < 



6'J l 26 



an iv r, u-vri o, o i o o 



0 7 



P c A L L , K , !< 0 , K 1 , K 2 

"REAOl 6, 1 ) A, t ,H, 0,RH ) , ! 

1 FORMAT! F6. A , A X , c r > . 2 , 5X , 3 ( n 6 . ° , S X ) ,"10.1 ) 

HOT If "(6,2) 

2 FORMAT! 1 H 1 / / i 6 X , 36 HI A 7 UP A L '' ' 1 D E S or CYIINOPR VIBRATION'/ 
112X.15HA 3 RAO I US' ( I N ) , 2 X , 1 S'U. = l c NGTII (IN), 

2 2 X , 1 P'HH = THICKNESS ( 1 N ) / 1 2 X , 1 0 HD = PPISSCN'S RATIO, 

— 3~2 X , 3 5 HR HO - OF N S 1 Tv (I n, S / C U I N ) / 1 2 X , 

431HP = YOUNO'S vnnoi ..US (IRS/SO ! N ) / 1 ?X , 2 1 HF P C C = FRF CUT NCY (K 

5 1 2 X , 321 IN ”=' NO . T C I P C U 0 F n R F N T I A L. I AVES/ 

6 1 2 X , ? 7 H M = NC. OF AXIAl H AL P -W A V F S /? •' X , 1 0 1 ‘ I N'RUT CATA/1 

HR I Tf ( 6, 3 ) A , I , H, 0, ROC , <= 

3 FORMAT ( 12X, A) A = , FA .4 , 5 X , AMI. = , F S . 2 , S v , A HN = ,F5.3/ 

1 1 2X , AMD = ,r 5.2,SX,swpor) = T T 5. 3 , ‘iX , 4HF = ,FIO.I/) 

WR I T P ( 6 , A ) 

A FOR M A T ( 2 7 X , 1 ] HfOJ 1 p UT [) A T A / 1 3 X ,1 F N , 7 X , 1 H m , <3 X , ] H K , 7 X , 5 H S P C F F , 

16 X , A H F R F 0 / ) 

P I = 3 » 1 A I 5 92 8 
TOl. = 0 . 00000 1 
rn 3;. o = 2, l 
00 P.5 L 1 = 1, 10 
v = L l-l 
P M= M 
RN = N 

K =A 6-p I Y R M / l 

«* B 1= ( H'X Y 2 ) /( 1 2.YAYY?) 











in 
1 6 
1 6 



20 



B = R l 

Y 1 1 =K. +11.-0) *no** 2/2. 

Y 1 2 =- { 1 . + 0 )*K*RNj/2. 

Y 1 3 =T)*K"" 

Y23=-( RN i R* ( P f i Y 3 + ( 2 .-0 ) yryy 26 (Ml 

Y 2 2 = ( 1 . +0 ) «K**2/2.+ 1 N*<*2 + n* { RJ* *2+2. *( I .-0) *K* *2 ) 

Y 3 3 - 1 . + R ■' { !< Y * ? + R N) v Y 2 ) Y Y ? 

K 0 = — Y 1 1 yy2?y Y3°-2. yv 1 2 * Y 2 3* y l 3 * Y 1 3 y yc> y y? 2 + y 1 2* y 2 * Y 3 3 <- Y2 3'*'*? y y 
K 1 = V 1 1 Y Y 3 3 : Y I } YV22+Y ??YY3 3-Y ] 3YY2-Y1 2 yy p-Y?? YY? 

K 2= - Y 1 1- Y22-Y33 
XO = --P. 2 5 
00 16 1=1,102 

f 1- v o Y Y 3 - !< 2 Y X p Y Y 2 4X1 Y X 0 - K 0 
F 2 “ 3 . Y X 0 Y Y ? - 2 . < •( ? Y V M K 1 
X ] - XP- r ' 1 /F2 

I F ( At SIX i-XO )-TO|. ) 5 3 , 1 6, 1 0 
XO=X 1 
CONT I OOF 

FRF 0 = Sn.’T ( ' 2 S ( xp ) Y F y 3«‘>.02/( ?HOY ( A + H ) yary ( ] ,-p^. 2 ) ) ) / ( 200°. ^-P 
y 2 = ■> B S ( X 0 ) 

X 3= SOR 1 ( v 2 ) 

WR I T F { 6 , 2 ! ) K V , R ‘ , K , X 3 , F R F 0 



A) 1 V G l. FVil 

7 l 
2 0 
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o, MOO 0 y A I K D AT C = or. 126 07/ 

f o R M A r ( 1 2 x , 2 ( r- 0 . y , AT ) , 2 ( F 0 . A , A X ) , F 7 . A ) 

CONT 1 MUF 
END 



. 
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NATURAL ««n = S OF CY! INOE y VIRRATin\| 

A - VADIMS (l\) 1. = I (NOTH ( I A' ) H - THICKNESS 

n = P O f 5 S n N ' S RUiG RHH = TENSITY (I.RS/CU IN) 

I = YOUNG'S MODULUS (LRS/SQ ID 
FREQ = FREQUENCY (KG) 

N = NO. Or C I p C lJ y F f-P C N T I A L WAVES 

Y = no. nr axial half- v.av'S 

"INPUT PAT A " 

A = 3. I 72 S' L = GO . OP H = 0.2 30 

n = 0.290 RHO = O.^T F = 29600000.0 



OUTPUT PAT A 



M 


T 1 


K 


SR 06 L 


FREQ 




JL 


A 


Zv Vc - 




2 . 0 


0, 0 


0.0 


0.06 83 


0.66 10 


2.0 


T . 0 


0.1661 


0. 07 98 


0.7720 


2. 0 


2.0 


0 . '<3 22 


C. 1 0 « 9 


1.0549 


2.0 


3.0 


0.4983 


0. 14 39 


1 . 4 4 C 7 


' 2.0 


A . 0 


0 . 6 6 4 4 


C. I9 60” 


l . 89 7 1 


2.0 


5.0 


0.8306 


0.24 80 


2 . 4000 


2.0 


6 . O' 


O'. 9 9 6 7 


0. 30 2 7 


2 . 9354 ' 


2.0 


p. 0 


1 . I 62 8 


0 . 35 R 3 


3. 4677 


2.0 


8.0 


1.3289 


0.41 43 


4. 000 O' 


2.0 


0 . 0 


1 . A 950 


0.4665 


4.5143 


3.0 


P.O 


0.0 


0 . n 3 ? 


1 .8693 


3.0 


! . 0 


0.1461 


0. 1058 


1 . 3950 


3.0' 


2.0 


0. 3322 


0. 20 38 


1 . 9 7 25 


3.0 


3.0 


0. 4983 


0. 2L 73 


2. 1025 


3.0 


4 . 0 


O'. 6 64 4 


0. 4 3 6 I 


2 .7346 


3.0 


3.0 


0. 8 7 06 


0. 2 ; , 00 


2.5157 


3.0 


6.0 


0.9967 


0. 2 3 89 


2.7906 


3. 0 


7.0 


1 . T 62 8 


0. 4 2 00 


3.1016 


3.0 


3 . 0 


I . <249 


0.45 5 6 


3.4412 


3.0 


9.0 


1 .4 950 


0. 39 2? 


3.3011 


4.0 


0.0 


0.9 


0 . 3 7 0 3 


3". 6 3 49 


4.0 


1 . 0 


0. 1 661 


0.4/16 


7 . 59 6 1 


4 .0 


? . 0 


9. '<3 22 


Oo 4 r 64 


3 . 6 3 2 1 


4.0 


3.0 


p.4964 


9.7’ 19 


3.6967 


4 . 0 


4 . n 


). 6 ' ■ r 4 


0 . 3 -> 11 


4.7347 


4 „ p 


6 . P 


A.SV • 


0 . 4 •* 4 1 


4. 5 p 1 2 


A A 

* e J 


• . o n 


0 , 9 9 7 


o.4i 0 r 


4 . p 4 5 


A * 0 


7 . 0 


1 . ] ' • 2 - 


0. 4 3 6ft 


4 . ' 1 7 6 


A oO 


0.0 


1 . 3285 


0. 45 6 4 


4 . M 6 6 


A * 0 


9. 0 


I .4960 


P . 4 7 0 5 


4 . 64 C 7 



( T N) 
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